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The Tabulation of Mathieu Functions 


Introduction 


The equation known as the Mathieu equation, which we write, following 
Ince, 


(1) d*y/dx* +- (a — 26 cos 2x)y = 0, 


and in which a and @ are constants, arises in applied mathematics in two 
main groups of problems. The first group, of which Mathieu’s problem of 
the transverse vibrations of a taut elliptic membrane is typical, consists of 
boundary problems relative to some partial differential equation, such as 
the wave equation, where the boundary is an ellipse—or possibly a hyper- 
bola. Here a fundamental requirement is usually that y shall admit the 
period 2x in x. With a prescribed value of 6, this restricts a to a discrete 
set of characteristic numbers, each associated, in the concrete problem, 
with the frequency of a normal mode of vibration or some analogous quan- 
tity. In the complete solution of this class of problem, solutions of the 
associated equation 


(2) _ @y/dx* — (a — 26 cosh 2x)y = 0, 


will usually also be needed. 

A typical example of the second group of problems is the modulation of 
a radio carrier wave. Here one is concerned with either a steady state or 
an initial value problem for an ordinary differential equation, and the 
independent variable x will frequently denote time. In this case both a 
and @ are prescribed. It is known that the general solution may be written 


(3) y = Ae*o(x) + Beo(— x), 


where A and B are arbitrary constants, ¢(x) is periodic in x admitting the 
period 27, and yu is a quantity depending on a and @ termed the ‘character- 
istic exponent’. Interest centres primarily upon whether yu is real or imagi- 
nary'—in the former case the solution is ‘unstable’, whereas in the latter it 
is ‘stable’, i.e. remains bounded for all (real) values of x. 

For the periodic solutions, 4 = 0, so that the curves in the a-@ plane 
corresponding to the characteristic numbers dissect this plane into regions 
such that adjacent regions correspond to stable and unstable solutions 
respectively. For both groups of problems, then, knowledge of the character- 
istic numbers as functions of @ is a primary need. 


The periodic solutions 


The periodic solutions of (1) which admit the period 2 fall into four 
classes, according as their Fourier series involve respectively cosines or 


1For real values of a and 6, » may be chosen, with no loss of generality, to be either 
purely real, or purely imaginary. 
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sines of even or odd multiples of x. They are defined by Ince as follows: 


(4.1) Cean(x, 0) = Lae cos 2rx, 

(4.2) 0) = sin(2r + 1)x, 
(4.3) &) = ABH? cos(2r + 1)x, 
(4.4) 0) = sin(2r + 2)x, 


corresponding to which the characteristic numbers are denoted in order by 
den, Denti, 


In every case satisfaction of the differential equation necessitates 
recurrence formulae connecting three (in the first, two) successive coeffi- 
cients. Elimination of the ratios of successive coefficients leads to an equation 
(involving a continued fraction) which expresses the relation between a 
and @. The continued fraction forms the basis of the technique developed 
by Ince and Goldstein for the computation of the characteristic numbers, 
and this technique seems generally to have been used by later computers. 

When the characteristic numbers are known, the ratios of the coefficients 
are calculable, but to make the coefficients definite, some normalisation 
rule must be chosen. 


Tables 


1. E. L. Ince, “Tables of the elliptic-cylinder functions,”’ R. So. Edinb., 
Proc., v. 52, 1932, p. 355-423. 


2. E. L. INcE, “Zeros and turning points of the elliptic-cylinder functions,” 
tbid., 1932, p. 424-433. 


1 is believed to be the only published table giving the numerical 
values of the functions themselves, and we therefore consider it first, before 
others chronologically earlier or later. 


The tabular contents of no. 1 are 
T. I, Characteristic numbers, p. 363. This gives doc1)5 and b1,1)6 for 
@ = [0(1)10(2)20(4)40; 7D]. 
T. II-XIII, p. 364-375. These give Fourier coefficients for céo.1)5(x, 6) and 
$€1(1)6(x, 8), to 7D, for the same values of @ as T. I. 


T. XIV-XXV, p. 376-423. These give the functions ceo 1)5(x, @) and 
$€i(1y6(x, 0) for x = 0(1°)90° and @=[1(1)10; 5D], along with second 
differences. 


An introduction, p. 355-362, defines the notation and the functions, 
outlines the methods of computation, and describes the contents of the 
tables. 
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The normalisation rule adopted is that, if y is any of the functions 
tabulated, then 


(5) f “ide 2. 


It follows that no Fourier coefficient exceeds unity in absolute value. The 
sign is chosen to secure that the values of the ce functions and the derivatives 
of the se functions shall be positive at x = 0. As a consequence of these 
rules it follows that, as @—>0, céo(x,0)—>1/,/2, while for m > 1, 
Cm(x, 0) —> cosmx and sé»(x, 6) — sin mx. 

Values of the functions for values of x in quadrants other than the first 
can be obtained from the tables, since the properties of the functions as 
regards symmetry and antisymmetry about x = 0 and x = + $7 are those 
of the trigonometrical functions to which they tend as @— 0. Values for 
the corresponding negative values of @ can also be obtained, since changing 
the sign of @ in equation (1) is equivalent to replacing x by ($x — x). It is 
to be noted that changing the sign of @ also, in effect, interchanges a,, and 
bm for odd (but not for even) values of m. 

Interpolation of the functions in the x-direction is readily accomplished 
by means of the differences provided, but interpolation as regards @ would 
seem not to be feasible at this interval, neither in the functions nor in the 
characteristic numbers and Fourier coefficients. 

T. 1 seems to have been completely overlooked in America; the present 
writer has encountered no reference to it in any American source. 

The tabular contents of no. 2 are 


T. I. Zeros, p. 431. This gives the zeros in degrees and decimals to 4D, of 
the functions cé2,1)5(x, 6) and sé3.1)6(x, 8), in the open interval 0 < x < 42, 
for 6 = 0(1)10(2)20(4)40. 


T. II. Turning points, p. 432-3. This gives the turning points of the functions 


C€1(1)8(X, 8) and Sée,1)6(x, 0), for the same range of x and values of @ as in 
Table I. 


The numerical results exhibit the asymptotic approach of ce», and Sé@m41 
as @ increases, which had been anticipated on theoretical grounds. The 
introduction contains formulae which give approximations to the zeros and 
turning points for large, and also for small, values of @. 


Mention must also be made of short tables of characteristic numbers in 


earlier exploratory papers by Ince, not all of which are incorporated in 
no. 1. The papers are 


3. E. L. INck, ‘‘Researches into the characteristic numbers of the Mathieu 
equation,” R. So. Edinb., Proc., v. 46, 1925, p. 20-29. 


4.—(second paper), ibid., v. 46, 1926, p. 316-322. 
5.—(third paper), ibid., v. 47, 1927, p. 294-301. 

In these papers Ince used as the Mathieu equation 
(6) d’y/dx* + (a + 16g cos 2x)y = 0, 


so that the don: and bens: of the later table no. 1 are interchanged. The 
sequence of characteristic numbers, in order of magnitude, is here 
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Qo0bybededsbs --- as against aob1@1bede2b3a3 --- in no. 1. It is to be inferred 
that Ince became convinced of the desirability of the changes of parameter 
from g to 0, and of sign, during the course of the computations for no. 1. 

In no. 3 Ince develops the continued fraction technique for the compu- 
tation of the characteristic numbers, exhibits it by numerical examples, 
and gives a table of numerical results. The table contains a» a; 6; b2 and a; 
(in the earlier notation) for g = [0(.1)1; 5D]. A diagram exhibiting the 
characteristic curves in the g-a plane is also given (p. 28). 

In no. 4 Ince explores more deeply an asymptotic formula for the 
characteristic numbers, valid when q is large. In a table (p. 321) doc1)5 and 
biqys (older notation) are given for g = [0(.5)5; 6D]. Since the continued 
fraction process is one of successive approximation, it is a great advantage 
to be able to calculate directly a good first approximation, and herein lies 
the great value of this asymptotic development. It appears, however, that 
this does not discriminate between a’s (or 6’s) adjacent in the (older) 
sequence, and is really an expansion for some intermediate characteristic 
number. 

In no. 5 this last point is further explored, and the fact that the asym- 
ptotic expansion relates to the characteristic numbers of half-odd-integral 
order is shown by numerical instances. The tables, (I-IV, p. 296) contain 
Qo 1, 61 Ay b2, G2 a3, and bs ag bs, for g = 2, 4.5, 8, 12.5, 18 (ie. 
(32g) = k = 8(4)24). 


Great as was Ince’s published contribution to the tabulation of the 
Mathieu functions, it represents a considerable curtailment of what he had 
planned, and indeed accomplished. He had worked out, to at least 12 
decimals, all the characteristic numbers which he published in no. 1, and 
also, to at least 10 decimals, a@o.1)4 and 61,1), for the values 64, 100 and 144 
of 6. His Fourier coefficients were originally calculated (using a normalisation 
rule different from that ultimately adopted) to 12 decimals, and there 
exists a manuscript table of the twelve functions as tabulated in no. 1, to7 
decimals, with second differences. The manuscripts have been loaned to 
the present writer by L. J. C. 


6. S. GoLpsTEIN, “‘Mathieu functions,’’ Cambridge Phil. So., Trans., v. 23, 
1927, p. 303-336. 


What seems to be the first systematic attempt of any magnitude to 
give numerical results for the periodic Mathieu functions was published on 
unpaged folding sheets at the end of this memoir. 

Goldstein writes the Mathieu equation in the form 


(7) @Py/dx* + (40 — 16g cos 2x)y = 0, 


so that he writes 4a where Ince writes a (and 4am instead of a, and 48, 
instead of b,,) and also uses the older parameter g. The tabular material 
consists of : 


T. I-V, which give for ceo, sé:, ce:, s€2 and ce respectively the values of 
ao, 81, a1, B2, and az, to 5D, along with the Fourier coefficients, for 
q = [.1(.1)1(.2)4, 5(5)30, 40(20)100(50)200; 5D]. 
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In the memoir Goldstein drew attention to the unsuitability of the 
older ‘normalisation’ rule whereby the coefficient of cos mx in Cém, or of 
sinmx in S€m, was taken as unity, and proposed the rule (5) above, except 
that he chose 2x as the value of the integral for the function céo, which 
then tends to unity as g tends to zero. Thus Goldstein's Fourier coefficients 
agree with Ince’s, except that for cep they are 2 times as great, while 
the characteristic numbers tabulated by Ince are four times those tabulated 
by Goldstein. Goldstein gives a very much wider range of q(@) than Ince, 
or indeed any other known table. The only values which overlap with 
Ince’s are 


Goldstein g=.5 1 2 3 4 5 
Ince 8 16 3 


Goldstein makes considerable use of the asymptotic expansion due to 
Ince, and in the memoir gives many more terms of it than Ince had calcu- 
lated. The asymptotic expansion is further investigated in 


7. S. GotpsteIn, ‘On the asymptotic expansion of the characteristic 
numbers of the Mathieu equation,’’ R. So. Edinb., Proc., v. 49, 1929, p. 
210-223. 


Here, after lengthy and complicated analysis, the dominant term of the 
difference between a, and an,; and between d@n4; and bn;; is determined. 
A second term is determined empirically and in T. III, p. 222, the values 
of b; — ao, b2 — ay, bs — a2, and by — a3 are compared with the values 
given by the two terms, for g = [.5(.5)5; 6D]. T. IV, p. 222, shows that 
these results are good for small values of g by comparing 6; — ao and Bz — a 
with the two terms of the asymptotic expansion, for g = [.1(.1)1(.2)2.6; 5D]. 


8. H. P. MuLHOLLAND & S. GopsTEIN, “The characteristic numbers of 
the Mathieu equation with purely imaginary parameter,” Phil. Mag., s. 7, 
v. 8, 1929, p. 834-840. 


In this paper the authors put g = is in (7) and give the following Tables: 
T. Ia, p. 839. ao and ae for s = [.02(.02).18; 6D]. 
T. Ib, p. 839. ao, 81, a3 and ® for s = .2(.2)2; and B: for s = [.2(.2).8; 6D]. 


It was found that a» and as, which for small s are real and distinct, 
become equal when s = .1836---, and for greater values of s (at least as 
far as s = 2) are conjugate complex quantities. Similar remarks apply to 
B: and ;. These facts explain the arrangement of, and what at first sight 
appear to be gaps in, the tables. 


9. K. Hipaka, “Tables for computing the Mathieu functions of odd order 
@), ces(x, 0), ses(x, 0), sé7(x, 0), and ce;(x, 6), and their derivatives,” 
Imperial Marine Observatory, Kobe, Japan, Memoirs, v. 6, no. 2, 1936, 
p. 137-157. 


The author uses (essentially) Ince’s notation and method of normalisa- 
tion. The Tables (p. 142-157) give biceyz and @1(2)7 for 6 = [0(.1)2.3; 7D}, 
along with the Fourier coefficients of the functions, B2,;; and A241, and 
also those of the derivatives, — (27 + 1)Beor4: and (27 + 1)Aor4:, to 7D. 
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The restriction to functions of odd order is apparently dictated by the 
requirements of a tidal problem for which the tables were computed. 
Hidaka asserts that the interpolation for @ is possible. A few tests by the 
present writer would indicate that while second differences may be adequate 
for the Fourier coefficients, the greatest possible accuracy in the character- 
istic numbers requires the use of fourth differences. 


10. J. A. Stratton, P. M. Morsg, L. J. Cuu, & R. A. Hutner, Elliptic 
. Cylinder and Spheroidal Wave Functions, New York, Wiley, 1942, xii, 127 p. 


This consists of a photographic reproduction of a paper by Chu & 
Stratton, with the same title, published in J. Math. Phys., v. 20, 1941, p. 
259-309, repaged as p. 1-51 in the book, to which have been added some nu- 
merical tables and an explanation of them, photographically reproduced 
from typescript. It was reviewed by H.B. in MTAC, p. 157-160. 

As far as the elliptic cylinder functions (Mathieu functions) are con- 
cerned, the equation used is 


(8) PS/dg? + (b — ccos*¢)S = 0, 


so that the relations between the parameters used here and those of Ince are 
(9.1) 6= 7/4, c = 

(9.2) a=b-c/2, b=a+ 28. 

The periodic solutions of (8) are defined by 

(10.1) Sei(c, cose) = cos ng, 

(10.2) Soi(c, cos¢) = Fy sin ng, 


where the primc indicates summation over positive integral values of of 
the same parity as /. It will be seen that Se and So correspond to ce and se 
in Ince’s notation. The normalisation rules adopted are that when ¢ = 0, 
Se and the ¢-derivative of So are each unity, i.e. 


(11) = 1, L'nF, = 1. 


It follows that the relations between these Fourier coefficients, D and F, 
and those of Ince, A and B, are 


(12.1) AQ = = 
(12.2) Bn = = By. 


The values of b corresponding to the solutions Se and So are denoted 
by } and Db’ respectively, and the correspondence between these and the 
characteristic numbers of Ince will be clear from the above. 

The Tables devoted to elliptic cylinder functions contain (p. 77-85) 
Docya and b’1(1)4 [5S], and the corresponding Fourier coefficients D or F 
[4 or 5D], for c = 0(.2)4.4, the values .1, .5(1)4.5 being intercalated. The 
greatest Fourier coefficient is given to 5S, and the others to the same 
number of decimals. 
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11. S. Lusxin & J. J. Stoker, “Stability of columns and strings under 
periodically varying forces,” Quarterly Appl. Math., v. 1, 1943, p. 215-236. 


Since the interest of these authors is in stability, characteristic numbers 
serve their purpose, and the tables at the end of the paper are restricted to 
these. The differential equation is taken in the form 


(13) @f/d? + (a + B cos 6)f = 0, 
so that the relations between their parameters and those of Ince are 


a = a/4, a = 4a 
B= —0/2(=—49), — 2. 


Corresponding to the sequence 
ao, bi, a, bo, a2, bs, 


of the characteristic numbers as used by Ince, Lubkin & Stoker employ 
the notation 


a(Co), a(C,), a(S), a(S2), a(C2), a(Cs), 
The Tables (p. 232-235) contain 


a(Co) --- a(S2), for 6 = [0(.2)4(.4)10(1)16(2)20; 5D], 
a(C2) --- a(Se), for = [0(-2)4(.4)6(2)20; 5D]. 


These Tables reproduce—to 5 instead of 7 decimals,—12 out of Ince’s 20 
values, and for the same orders. Goldstein’s results have been used for 
comparison, where the overlap occurs (but see ‘Discrepancies’, below). 

It is known that there cannot co-exist, for the same values of a and @, 
two independent solutions of (1) both admitting the period 2x. In other 
words, the second solution corresponding to cén(x,0) or Sém(x,@) is not 
periodic. This second solution is rarely of physical importance, but is 
mathematically interesting—and elusive. One attempt to define and provide 
a basis for calculating a second solution is known. 


12. S. GoLpsTEIN, ‘‘The second solution of Mathieu's differential equation,” 
Cambridge Phil. So., Proc., v. 24, 1928, p. 223-230. 


The table (p. 230) gives coefficients b{";{” for r = 0(1)8 and m = 0(1)8 
for the second solution corresponding to ce:(x,qg) with g = 1 [8D], the 
definition being 


(14) iny(x, g) = 2xcex(x, g) + sin (27 + 1)x. 


In this, the coefficient b{7{{” is the contribution to the total coefficient of 
sin(2r + 1)x due to A%,4:. 


13. Puttip M. Morse & Peart J. RUBENSTEIN, Tables of Elliptic Cylinder 
Functions. Calculated in 1940, 3 p. of formulae and 17 p. of tables. Hecto- 
graphed from typescript. 

The writer is indebted to R.C.A. for a photoprint of these tables, which 
has enabled the following account to be given. The tables form part of 
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the basis for the corresponding tables in no. 10, which have in some respects 
superseded them. 

In connection with the periodic solutions Sen(c, cos¢) and Son(c, cos¢), 
the tables give, in addition to the separation constants 6, and b’», the 
quantities N,, and N’,, defined by 


‘2a 
(15.1) Stine f = 
(15.2) N'n = f = 


where the double prime indicates summation over all positive integers and 
zero, of the same parity as m, together with the factor 2 for the terms in 
which n = 0. ; 

These tables are notable, however, for the fact that they give some 
numerical information as regards soiutions of the associated ‘hyperbolic’ 
equation. The differential equation, equivalent to (2), is taken in the form 


(16) @PRidy? + (c cosh*y — b)R = 0. 
Clearly one solution of (16), with b = bn, is 
(17) Sén(c, cosh y) = cosh ny. 


This solution converges slowly if y is not small. There is, however, an 
alternative solution of (16) in the form of a series of Bessel functions, 
which has convenient asymptotic properties, and it is remarkable that the 
coefficients in this expansion are multiples of D7. In fact, if 


(18) Rem(c, 2) = J, (cz), 
then 
(19) Sém(c, cosh = V(22)AmRem(c, cosh 


The satisfaction of the differential equation (16) by the series in (18)— 
with z = cosh Y—depends upon the differential equation and recurrence 
relations satisfied by J,(cz). These are shared by the functions of the second 
kind, N,(cz) in the notation of JAHNKE & Empk, or Y,(cz) in the notation 
of Watson and the BAASMTC. Consequently there is an independent 
second solution of (16) 


(20) Reém(c, cosh = cosh y). 
With y = 0 in this, u, is defined by 

(21) Ren(c, 1) = um; 

and also C,, and ym by 

(22) Cn = + COtYm = 


Corresponding to Som, with the separation constant b’,, the functions Ron 
and Rom, and the constants Y'm, are similarly defined. 
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The first group of tables gives },, to 5S, Nn, to 5S, vV(2x)Am, to 5S, 
V(27) um, to 4S, C,,, to 4S, ym(°*), and the Fourier coefficients D7, to 4D, 
for m = 0(1)4, c = O(.1).2(.2)4.4 and 4.5, and the corresponding values 
for the odd functions for m = 1(1)3 and the same values of c. 

The second group of tables is easier to describe. It consists of values of 
Seocya(c, cos and Soi14(c, cos ¢) for c = 0(.5)4.5 and ¢ = 0°(5°)90°; 
[4D]. 

A quite casual examination has shown numerous discrepancies between 
these tables and those of no. 10, of which a few are not small. The reliability 
is therefore not high, but they contain the results of considerable numerical 
exploration, provide results that might be adequate for some purposes, 
and could be very valuable in giving first approximations to anyone who 
wished to achieve greater accuracy. 


14. W. G. BickLey, Characteristic Numbers of the Mathieu Functions. 
Unpublished manuscript. 


The present writer has for some years gradually compiled a table of 
values of the characteristic numbers for the values of k = ¥@ = 0(1)10. 
The calculations have been carried to sufficient figures to determine twelve 
decimals of the characteristic numbers. Results are complete as far as a; 
while some for bs have been obtained. 

Some graphs showing the dependence of the Fourier coefficients BE";{” 
upon g are given in 


15. I. Lotz, “‘Korrektur des Abwindes in Windkandlen mit kreisrunden 
oder elliptischen Querschnitten,” Luftfahrtforschung, v. 12, 1935, p. 250-264. 


Lotz follows Goldstein no. 6 and gives, in figs. 17-19, p. 259, graphs 
of the first four or five Fourier coefficients for se:, ses and ses, for the range 
of 0 <q < 20. 


Although no numerical tables are given (readers being referred to no. 6 
and no. 8 above) reference must be made to 


16. E. JAHNKE & F. Empe, Funktionentafeln mit Formeln und Kurven, third 
ed., Berlin, 1938; or American reprint, New York, Dover Publications, 1944. 


P. 283-295 give a synopsis of the formulae, following the notation of 
Goldstein no. 6. 

Fig. 166, p. 286, gives the a — g curves for the ranges — 13 < g < 13, 
— 6 < a < 14, for the characteristic numbers from ap to 67. 

Figs. 167-177, p. 288-293, give various graphs of the five functions, 
ceo(x, g) to see(x, g) while figs. 178-181, p. 294, 295, exhibit graphically the 
results of Mulholland and Goldstein no. 8. 


Non-periodic solutions of the Mathieu equation depend on the char- 
acteristic exponent—cf. equation (3) supra. Some values associated with 
this exponent, and also some graphs of non-periodic solutions, are given in 


17. J. G. Bratnerp & C. N. WeyGanpt, “Solutions of Mathieu’s equation 
—I,” Phil. Mag.,s. 7, v. 30, 1940, p. 458-477. 
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The Mathieu equation is used in the form 
(23) d’y/d? + «(1 + kcost)y = 0. 


Two fundamental solutions, g(t), h(t), are defined, such that g(0) = 1, 
g’(0) = 0; (0) = 0, h’(0) = 1. It is shown that 


(24) g(2r) = h'(2r) = cos = BD. 


Numerical values of g(#) and A(#) were determined for a range of values 
of e and k by numerical integration of the differential equation, and hence } 
was determined. Sample graphs show g(t), A(#), and their derivatives, for 
e=9,k = .9. In 


T. I, p. 469, values of 5 are given for k = .1(.1).6, .8(.1)1 and 3v « = [0(1)6, 
8, 9; 3D]. 


Notation 


The conflicting notations used by the various writers are confusing. 
After considerable thought and discussion, the notation for the canonical 
form of the Mathieu equation to be used in a forthcoming account of the 
functions and their applications is 


d*y/dz? + (a — 2q cos 2z)y = 0 


with Ince’s notation for the characteristic numbers and for the functions, 
and with Ince’s normalisation rule. The reason for the introduction of the 
coefficient 16 by Mathieu has lost its force, and the mixture of Greek and 
Roman letters in (1) cannot easily be defended. Physically, the square 
root of the numerical value of g is significant, and this is also useful in 
defining solutions of the associated hyperbolic equations. It is proposed 
to use k such that k? = |g]. 

Further details, and definitions of a set of functions satisfying the 
Mathieu and its associated equations, will be given elsewhere. 


Discrepancies 


As far as the writer is aware, no error has so far been reported in no. 1 or no. 2. One 
error in no. 6 has been reported, by Goldstein himself (in no. 7). The value of a; for g = 2, 
should be 2.33382 and not 2.33817. 

In the overlapping cases (c = 2, 4, @ = 1, 4) the separation constants given by no. 10 
are in agreement with the characteristic numbers of Ince. The labour of comparing the 
Fourier coefficients (differently normalised) has not been undertaken. 

The overlapping values in no. 11 have been compared with those of no. 6 and no. 7. 
As regards no. 6 the following discrepancies occur: 


Goldstein no. 6 Lubkin & Stoker no. 11 
g=0.3 ao = — 0.50534 B=1.2 a(Co) = — 0.50535 
q=04 a = — 0.77897 B=1.6 a(Co) = — 0.77898 
g=0.2 a= 0.55406 B=0.8 a(Si) = 0.55906. 


Since Lubkin & Stoker had compared their results with Goldstein’s, one may presume 
that they have checked the end figures of the first two of these. Differences indicate that 
in the third of the above, Goldstein is correct. 


Cc 


, 6 

6 

6 

6 

6 

6 

6 

6 

In th 

error, 

value: 

evide 
Cc 
In th 

25.02 

B 

Addit 
F 
| 
nik, ( 
E 

: Press 
I 
Univ 

Impe 

LON 

\ 
I 

Wa’ 


1e 
d 


10 
he 


ne 


RECENT MATHEMATICAL TABLES 


Comparison of no. 1 with no. 11 yields the following: 

Ince no. 1 Lubkin & Stoker no. 11 
6=40 az = — 20.20794 08 a(C2) = — 5.05198 (= 20.20792 + 4) 
2 &= 9.14062 77 a(C;) = 2.28515(= 9.14060 + 4) 
4  10.6710271 a(S3) = 2.66777 (= 10.67108 + 4) 
bh = 13.55279 65 a(S;) = 3.38817 (= 13.55268 + 4) 
@=12 22.9721275 a(C,) = 5.74803 (error) 
6=32 bs = 26.10835 26 a(Cs) = 6.52721 (= 26.10884 + 4) 
6=40 bs = 22.3321485 a(Cs) = 5.58302 (= 22.33208 + 4) 
6=40 a;= 41.34975 44 a(Ss) = 10.33749 (= 41.34996 + 4) 
6=40 be = 41.43300 52 a(Ss) = 10.35813 (= 41.43252 + 4) 


In the above, differences indicate that Lubkin & Stoker’s value of a(C,) for 8 = 6 is in 
error, and should be 5.74303. In the other cases the discrepancies are too small for the 
available differences to discriminate with certainty, but in view of the fact that all Ince’s 
values were worked out to 12 decimals by a process of successive approximation, and bear 
evidence of careful checking, they are prima facie the more reliable. 

Comparison of Hidaka no. 9 with Ince no. 1 shows two end-figure discrepancies 


Ince no. 1 Hidaka no. 9 
@=2, B® = 0.12413 61 0.12413 60 
@=1, as = 25,02085 43 25.02085 44. 


In the case of the second of these, the present writer agrees with Ince, having obtained 
25.02085 43454 5... 


Comparison of Hidaka no. 9 with Goldstein no. 6 shows one discrepancy. 
Goldstein no. 6 Hidaka no. 9 
g=0.2 — BY = 0.16171 6=1.6 — BY = 0.16171 80 


Comparison of Hidaka no. 9 with Lubkin & Stoker no. 11 confirms that a(5S;) for 
8B = 0.8 should be 0.55406. It yields also the following. 


Hidaka no. 9 Lubkin & Stoker no. 11 
@=1.2 bs = 9.06485 47 B=0.6 a(Cs) = 2.26622 (= 9.06488 + 4) 
@=2.0 bs = 9.14062 77 B=1.0 a(Cs) = 2.28515 (= 9.14060 + 4) 


Additional References 


For collected accounts of ranges of the theory of solutions of the Mathieu equation, 
reference may be made (in addition to nos. 6, 10, and 16 above) to 

M. J. O. Strutt, Lamésche, Mathieusche, und verwandte Funktionen in Physik undTech- 
nik, (Ergebnisse der Mathematik, v. 1, part 3), Berlin, 1932. 

E. T. WaittaKer and G. N. Watson, Modern Analysis, Cambridge, University 
Press, fourth ed. 1927; Amer. reprint, 1943. 

E. G. C. Poote, Introduction to the Theory of Linear Differential Equations, Oxford, 
University Press, 1936. 

W. G. BICKLEY 

Imperial College of Science and Technology, 
LONDON, S. W. 7, England. 
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195[A, D].—J. G. BEcKERLEY, “The calculation of arg T (ta + 1),”” Indian 
J. Physics, v. 15, and Indian Assoc. for the Cultivation of Science, Proc., 
v. 24, 1941, p. 229-232. 16.5 X 23 cm. 


T'(z) for complex values of z has been tabulated to a rather limited extent. The tables of 
WALTER MEISSNER, 1939, -were reviewed in MTAC, p. 177; in his Tables of the Higher 
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Mathematical Functions, v. 1, Bloomington, Indiana, 1933, p. 269f., H. T. Davis tabulated 
1/T (re), to 12D, for r = — 1(.1) + 1, @ = 0°, 30°, 45°, 60°, 90°, 120°, 135°, 150°; in 
C. P. Weis & R. D. Spence, J. Math. Phys., M.1.T., v. 24, Feb. 1945, p. 61, there is a 
table of I'[(3 + ia)/4] giving the argument and modulus to 4D, a = 1(1)5; JAHNKE & 
Empe’s reference, p. 21 of the 1943 edition, to a certain GINZEL table of (x + 7y)!, is without 
foundation in fact. 

Beckerley gives a 4-place table of arg T'(ia + 1) fora = [0(.1)2; 4D], with indications of 
the methods of calculation. This function occurs frequently in problems involving positive 
energy hydrogen functions.! The author expresses the hope that ‘‘this table will be of use 
to physicists who are engaged in numerical calculations involving the ‘Coulomb phase 
factors’ which occur in the continuous spectrum wave functions of hydrogenic atoms.” 

x A. 


1A, SOMMERFELD, Wave-Mechanics, London, 1930, = G. Gamow, Structure of 

Atomic Nuclei and Nuclear Transformations, Oxford, 193 psy IX, uation 15, p. 163; 

oat ~~; a Fundamental Principles of Quantum Mechanics, New York and London; 
Pp 


196{A, F].—I,. J. KavAn, Rozklad VSetkjch Ctsel Celjch od 2 do 256 000 v 
Prvocinitele, and on the second title-page: Tabula Omnibus a 2 usque ad 
256 000 Numeris Integris Omnes Divisores Primos Praebens (Observa- 
torium Publicum, Star4 Dala, Czechoslovakia), Prague, 1934, xi, 514 p. 
28.3 X 29.4cm. See RMT 71, Scripta Mathematica, v. 4, 1936, p. 338. 
I,. Factor Tables giving the Complete Decomposition into Prime Factors of 
All Numbers up to 256 000..., with prefaces, English by B. STERNBERK, 
and Latin by K. Petr, and an introduction by ArTHUR BEER, 
London, Macmillan, 1937. xii, 514 p. 28.3 X 29.4 cm. 42 shillings. 
II. BAASMTC, Mathematical Tables, volume V. Factor Table giving the 
Complete Decomposition of all Numbers less than 100,000 prepared inde- 
pendently by J. Peters, A. LopGE and E. TERNouTH, E. GIFFORD . . ., 
London, 1935, xv, 291 p. 21.6 X 27.9 cm. Published at 20 shillings; out 
of print. For brevity these tables will be referred to in what follows as 
Table I and Table II. 


Many mathematical problems arising out of the war effort are concerned with the de- 
velopment of new expansions for very special functions. These expansions, in most cases, 
have rational coefficients. To maintain accuracy in the combination and checking of these 
expansions it is often best to keep denominators in factored form and, to reduce the rational 
numbers to their lowest terms, it is convenient to know the factors of the numerators also. 
In short, a good computer need not be a number theorist to profit by using a fair-sized 
factor table at appropriate times. 

A number of small factor tables will indeed be found in machinist’s handbooks and small 
volumes of collected tables. These usually extend to about 10000 and often omit multiples of 
2, 3, and 5, the very numbers most frequently met with in work of the sort described above. 
These embarrassing omissions are a practical necessity in the great factor tables of the first 
ten millions, which are intended after all for number theorists. The two tables under review 
are unique in giving all the prime factors of every number in their respective ranges and at 
the same time being sufficiently extensive to take care of the 5- or 6-figure work most fre- 
quently required. A comparison of the two tables may be made as follows. Table I is of 
course more than two and a half times as extensive as Table II. However, this does not 
mean that the former is twice as useful as the latter. Anyone who inspects a much used 
copy of a factor table cannot fail to notice the frayed and soiled edges of the early pages and 
the comparatively new appearance of its later pages. 

The two tables are well arranged, though quite differently. Table I is arranged like 
an ordinary table of logarithms, consecutive numbers being in adjacent columns which are 
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headed 0, 1, ---, 9, with 50 lines to the page. The lines are numbered at both ends to facili- 
tate entering the rather wide (23 cm.) table. Table II is arranged so that consecutive numbers 
occupy the same column and adjacent lines, an opening being devoted (in most cases) to 
700 numbers, as compared with 1000 in Table I. The left page of the opening gives the fac- 
tors of numbers whose last two digits are 00, 01, ---, 49. The average computer will no doubt 
find the arrangement of Table I more familiar and natural. When the factors of a long series 
of consecutive numbers are needed (as in some number theory problems) Table II is more 
convenient to use. In case the number to be factored is actually a prime both tables print the 
number in bold-face type, thus assuring the user that he has entered the table correctly. 
This is a decided improvement over earlier tables in which primes are represented by a dash. 

The printing of Table I is decidedly inferior to that of Table II. There is much irregu- 
larity in type and inking. The exponents of the primes are unnecessarily small and in some 
cases (e.g. N = 123725) are of two sizes. The printing job in Table II is beautifully done, 
the primes being in old type and the exponents in new type. 

Both tables are extremely reliable, especially Table II, which is based on three inde- 
pendent manuscripts and an intense program of proof reading and comparison with previous 
tables. The history of factor tables shows a preponderance of printer’s errors over author's 
errors. Table I, which took 17 years to prepare, was verified in manuscript by actually 
multiplying together the decompositions given. Had this been done on the proof sheets, the 
more serious of the two following errors might have been avoided: 


p. 32 N = 15280, for 24-3-191, read 24-5-191; 
p. 39 argument left column, for 1800, read 1870. 


These errors were discovered by J. C. P. MILLER. 

For the sort of work encountered in practical problems the numbers to be factored 
seldom contain large prime factors. Such numbers constitute a quite small minority, so 
that a small factor table devoted exclusively to such numbers may nevertheless have a 
considerable upper limit. Such a table has been published by Cunningham.! It gives the 
factorization of all numbers up to 100000 having no prime factor in excess of 11 on nine 
small pages, 1196 entries in all. 

D. L. 
Note sy S. A. Jorre: In the Czech introduction by KavAn there are two errors in quota- 
tions from CayLey, Collected Math. Papers, v. 9% p. 462-463. On p. V, box, factors of 391 
should be 17-23, mot 17-33; Cayley was correct and Kavan — the factor 23 incorrectly. 


On p. VI, first box, the number corresponding to 297 should be 180, mot 198; here Kavan 
copied Cayley’ s error, p. 463 


_ 1A. J. C. CUNNINGHAM, Quadratic and Linear Tables, London, Hodgson, 1927, p. 162- 


197[C, D].— Natural and Logarithmic Haversines, ‘‘arranged from Bowditch, 
American Practical Navigator, 1938 Edition, Table 34. By permission 
of the HyproGRAPHIC OFFICE, United States Navy Department.”’ New 
York, Macmillan, 1943, 38 p. 14.2 X 21.2 cm. Paper cover, 30 cents. 


This publication is misleading in more than one respect. In the source indicated, in the 
above quoted footnote, the five-piace table is of logarithmic and natural haversines, 0° to 
120°, at interval 15’’; 120° to 135° at interval 30”; 135° to 180° at interval 1’. The natural 
haversines are in black-face type throughout. What is “‘arranged”’ is simply to abridge this 
table to every minute of arc throughout, five degrees to the page, with constant uniformity 
in face of type. Hence this does not justify the title-page legend, Natural and Logarithmic 
Haversines by Paut R. RIDER . . . and CHARLES A. HUTCHINSON, even though this table 
also occupies p. 187-222 of Navigational Trigonometry, published by these authors in 1943. 

The other criticism is of an error made by the Hydrographic Office, rather than by 
these authors. Nathaniel Bowditch (1773-1838) never had anything to do with such tables 
as those under review. In 1844, Nathaniel’s son, J. I. Bowditch (1806-1889) edited a volume 
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of tables taken from the New American Practical Navigator of his father. These were called 
Bowditch’s Useful Tables (247 small p.) ; of this volume at least a score of editions had been 
published up to 1932. The editions after 1868 were published by a department of the Govern- 
ment, now called the Hydrographic Office, and after 1880 the word New was dropped from 
the earlier title New American Practical Navigator. Up to 1903 neither the Useful Tables 
(“Bowditch’s’”’ had been dropped from the title), nor the American Practical Navigator, 
contained haversine tables. But both of these works had such tables in their 1911 edition. 
Furthermore, the Hydrographic Office in that year published the volume, with the partial 
title-page, ‘Useful Tables from the American Practical Navigator by Nathaniel Bowditch” 
(427 large p.), even though no Bowditch had anything whatever to do with most of the 
contents of the volume, including the haversine tables. 

The most elaborate published table of haversines is that of J. C. HANNYNGTON, Haver- 
sines, Natural and Logarithmic, used in Computing Lunar Distances for the Nautical Almanac, 
London 1876, 327 folio p. This is a 7-place table (except for the first 36° of the logarithms), 
0 to 180°, log haversines at interval 15”, and natural haversines at interval 10’. For the 
range 0 to 125° or 135° a similar excellent table was published earlier, RICHARD FARLEY, 
Natural Versed Sines from 0 to 125°, and Logarithmic Versed Sines from 0 to 135° ---, 
London, 1856, 90 p. And yet earlier another 7-place table of haversines, at interval 10”, 
from 0 to 120°, James ANDREW, Astronomical and Nautical Tables with Precepts for finding 
the Latitude and Longitude of Places ---, London, 1805, T. XIII, p. 29-148. There is a 
copy of Andrew’s very rare work in the New York Public Library. 

A. 


198[C, E].—J. R. Hutt & R. A. HULL, “Tables of thermodynamic functions 
of paramagnetic substances and harmonic oscillators,’ J. Chem. Phys., 

v. 9, 1941, p. 465-469. 19.5 XK 26 cm. 

This paper, which was not mentioned in MTAC, p. 119, contains tables of the Einstein 
functions 2/Z, In Z, 2*e*/Z* where Z = e* — 1. The tables give 5S for values of z between 
0 and 14.4 at intervals ranging from .02 to .8. There is a misprint in the heading of the fifth 
column of p. 468. Tables are also given for Q and In Q, where 


i 
Q= 2 e; 
m=—j 
j = 3 (1) 2} is the angular momentum quantum number for the substances considered, and x 
ranges from 0 to 6. Actually x = [0(.02).2(.04).28(.02).32(.04).4(.1)4(.2)6; 5S]. 
H. B. 


199[E].—P. B. Wricut, ‘Resistive attenuator, pad and network, theory and 
design,” part 3 of a 4-part paper, Communications, v. 25, Jan. 1945, p. 57, 
58, 60. 19.5 X 27.1 cm. Compare RMT 174, p. 358. 

Tables for 20 log e(e% = k® > 1) = 0(.01).2(.05).4(.1)4(.5)30(1)60(5)140, 150, of (a) 
sinh 20 = (k* — 1)/2k?; (b) tanh? @ = [(k? — 1)/(k* + 1)]*; (©) (1 — e*) = — 1)/k; 
sinh? @ = (k* — 1)?/4k; (e)' cosh 20 = (k* + 1)/2k*; (f) tanh? 30 = [(k — 1)/(k + 1)]*; 
(g) csch 20; (h) coth? 6; (i) 1/[1 — e~]; (j) csch* @; (k) sech 26; (1) coth? $0. These tables are 
to 5-9S. In part 1, of this 4-part paper, Aug. 1944, v. 24, no. 8, p. 52, 54, 56 there are tables 
for the same range of argument for (a) In k, 3 In k, 2 In k; (b) k, 1/k, k*, 1/k*; (c) (k — 1); 
2(k — 1), (& — 1)*; (d) — 1), — 1). 

1 The right-hand member is given incorrectly on each of the three pages as (k? 4- 1)/2k?. 


200[E, M].—F. StAsiemn & R. ScHLAFER, ‘“‘Numerische Berechnung von 
y(x) --f e“dt,”” Z. angew. Math., v. 23, Feb. 1943, p. 59-61. 
21.5 X 27.7 cm. 
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Here is a table of y = e* f” edt, for x = (0(.1)10; 4D], which satisfies the differen- 


tial equation y’ + 2xy = 1. There are no references to other tabulations of this function, 
such as by TERAZAWA, and MILLER & Gorpon, that of the latter being considerably more 
extensive than the one under review; see MTAC, p. 323. 


201(L].—R. P. Batpwin, “Tables of functions used in determinations of 
stellar ionization temperatures,” Northwestern Univ., Dearborn Ob- 
servatory, Annals, v. 4, part 14, 1940, p. [3]. 

Here is a table of f” x%dx/(e — 1), mostly to 4 or SS, for xs = 1(.5)12.5. Compare 
0 
MTAC, p. 140, 46, and p. 189, RMT 152. 


202[L].—K. E. Bissuoprp, ‘Stress coefficients for rotating disks of conical 
profile,” Am. So. Mech. Engrs., J. Appl. Mech., v. 11, 1944, A-8—A-9. 
21.6 X 26.6 cm. 


In the calculation of stresses it is found to be advantageous to use both the fundamental 
solutions of the hypergeometric equation 
x(1 —x)P” + [ce — +0 + 1)x]P’ — abP = 0, 
a+b =1, ab =o —1, c =3 for the point x = 0; and also those for the point x = 1. 
In equation (15) the power series for t#F(a + 2, b + 2; 3; #) is given incorrectly but the 
calculations for the tables seem to have been made with the correct formula. The tables 
give the values of Pi, —P;', Ps, —P:2', for ¢ = 0.3 and 0.36, x = [0(.01)1; 6S] and it is 


thought that they are accurate to within 5 parts in 2,000,000. The functions P are defined 
as follows: . 


P(x) = Si(x) = F(a, b;3;x) = C2S2(¢) + CS2(t),t = 1 — x; 
P2(x) = S2(t) = #@F(a + 2,b + 2;3;t) = CsSi(x) + C8i(x); 
Si(x) = — + 1)b(b + 1)Si(x) Inx + — + + 


— — — 2) 4); 
n=2 
= Via + — 3) + + — 3) — — 2) — V(n) 
Salt) = 1 — abt — — 2)!] — + + 1)S2(0) In t; 
n=2 


Vv, = Via+n—1)+ — 1) — — 2) — 
C; = (2/x) Sin (ax) /[a(a + + 1)), 
C2 = (2/2) Sin (ar) [: vio) | — Cos (az), 


i 1 1 1 
= 


a—-2 a a+i 
y = Euler’s constant. ¥(x) = I'’(x + 1)/P'(x + 1). 


— 2y — 2¥(a) — x Cot (ar) | 


Calculations were actually made for ¢ = .24(.03).36. Copies of the unpublished tables 
may be had on application to Fairbanks, Morse & Co., Beloit, Wisconsin. 
The tables also include values of the stress coefficients p:, p2, ps, G1, 92, s- 
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203[L].—A. Biocu, On the temperature coefficient of air-cored self-induct- 
ances,” Phil. Mag., s. 7, v. 35, Oct. 1944, p. 704. 16.9 & 25.3 cm. 
Tables of Fi = (x2/8)(1 — 2), Gis = (x*/8)(1 — x = [0(.1)10; SD], where —@, 

= = or = 1 —2W(x)/ 

[xX (x)], 2 = 1 — 4Z(x)/[xV(x)]. X, V, W, Z, W/X, and Z/V have been tabulated by H. G. 

SavinceE, Phil. Mag., s. 6, v. 19, 1910, p. 55f.; compare MTAC, p. 256. 


204[L].—H. Bucuuotz, Ausbreitung der Schallwellen in einem Horn 
von der Gestalt eines Rotationsparaboloides bei Anregung durch eine 
im Brennpunkt befindliche punktférmige Schallquelle,” Amnnalen der 
Physik, s. 5. v. 42, Mar. 1943, p. 433-435. 12.9 X 20.5 cm. 


am® (it) /a¢ = (ig). For r < 0, ¢ > 0, or for r > 0,¢ <0, 

> 0, ¢ > 0, or forr < 0,¢ <0, 


my (ig) = { ( - +) (2V (78) 
=0 


where y(¢) is a given function of ¢ and X. 
Of the function (2/x)4m/< (i¢) there is a perspective graph 0 < ¢ <6, —3 <7 <+3; 
and a table ¢ = 1(1)6, r = [—3(.5) + 3; 5 or 6S]. Also there is a table of the first three 


zeros, 71’, 72’, 73’ of ms (it), ¢ = 0(1)6; 71’ to 5D; 72’ to 4D, 73’ to 2D. There are also graphs 
of the zeros < r'< +7,0<¢ <6. 


205[L].—H. Bucuuorz, ‘‘Die Ausstrahlung einer Hohlleiterwelle aus einem 
kreisférmigen Hohlrohr mit angesetztem ebenen Schirm,”’ Archiv. f. 
Elektrotechnik, v. 37, March 1943, p. 160, 161, 163. 19 & 27 cm. 


There is a table of J,(jo,.) for s = 1, 2, 3 and » = [.5(.5)20.5; 7D]; for s = 1 and 
vy > 10.5, s = 2 and » > 16, the values of the function are all zero. The graphs of J,(jo,s), 
s = 1, 2, 3, are for 0 < » < 13. On p. 163 is a table of the exact values of D,(p), for 
s = 0(1)7, p = 0(1)7, s + p > 7; and of C,(p), for s = 0(1)6, p = 1 (1)7, s + p succes- 
sively 1, 3, 5, ---, 13, where 


(2p + 1)(2p + 2) --+ (26 + A) 
(2p + 20 + 1)a! 
(0) = 


The expression for D,™ as given by Buchholz in (5.5b) is incorrect; for (2p + 1), read 
(2p + 1). 


R. C. A. 


206[L].—J. Cossar & A. Erp&ty1, Dictionary of Laplace Transforms. 
Admiralty Computing Service, Department of Scientific Research and 
Experiment, London; Part 1, no. SRE/ACS.53, 1944, 42 leaves; Part 2A, 
no. SRE/ACS.68, Dec. 1944, 49 leaves; Part 2B, no. SRE/ACS.71, 
Feb. 1945, 56 leaves. 20.2 X 33 cm. Mimeographed on one side of each 
leaf. These publications are available only to certain Government 
agencies and activities. 

Of this monograph, Part 1 contains an Introduction, and sections on: notes and ab- 
breviations, general formulae, short table of Laplace transforms, bibliography. Only the 
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unilateral Laplace transformation is considered and the Dictionary consists of pairs of func- 
tions f(t), @(p), connected by the relation 


ery wat = Lis: 


called the Laplace transform of f(t). In early days ¢(p) was called the generating function of 
f{( and later writers spoke of Abel’s generating functions. The extension of the term generat- 
ing function from series to integrals was a natural development of the theory of generating 
functions of Lagrange and Laplace. The functions f(¢), ¢(p) are classified according to f(t) 
in Part 2, and according to ¢(p) in Part 3 (which has not yet been published). In Part 2A 
¢(p) is given for each f(t) of successive groups of functions: rational, algebraic, powers with 
arbitrary index, jump- and step-, exponential, logarithmic, trigonometric, hyperbolic, and 
composite elementary functions. In Part 2B, f(¢) includes the following forms: Bessel 
functions, modified Bessel functions, products of Bessel functions, Bessel integral functions 
and Fresnel integrals; Kelvin and Struve functions; sine, cosine, exponential and logarithmic 
integrals; Legendre, Gegenbauer, Jacobi, Hermite, and Laguerre polynomials; parabolic 
cylinder functions, Bateman k-function, and Whittaker functions; Legendre functions, 
Gauss’s series, general hypergeometric series, functions of two or more variables; and theta 


. functions. In many cases $(p)is expressed by means of the various types of Bessel functions, 


confluent hypergeometric functions, the incomplete gamma function and other functions of 
modern analysis. 


No attempt is made at completeness. Not all available sources have been consulted ; 
and even from the works consulted not all formulae have been included. Formulae which 
need much explanation, or which involve rare functions, have been omitted. Sometimes 
only a few typical examples are selected from an extensive list of related formulae. Then 
reference is made to the list where more involved formulae are found. 

It is pointed out in the Bibliography that the most extensive publiched list of Laplace 
transforms in existence is in N. W. McLacHLan & P. HuMBERT, Formulaire pour le Calcul 
symbolique (Mémorial des Sci. Mathématiques, v. 100), Paris, 1941, 67 p. 

Each of these Dictionaries should be very useful. 

8. CA. 


207[L]—NYMTP, Jacobi Elliptic Functions, Washington, D. C. 1942 
34 hektographed sheets and six mimeographed sheets, printed on one 
side only 21.5 X 35.5cm. Not available for general distribution. 


These tables were described MTAC, p. 125-126. 


208[L].—NYMTP, Table in J. R. Whinnery & H. W. Jamieson, ‘‘Equiva- 
lent circuits for discontinuities in transmission lines,’’ Inst. Radio 
Engrs., Proc., v. 32, 1944, p. 114. 21.8 K 28 cm. 


The table is of Hahn’s function 
So(a) = > sin? (nwa) 
a = [.01(.01)1; 4—5S]. There are also aan 0 <a <1, of So(a), and of the associated 
function 
Sp(a) = p*sint — p = 1(1)8. 
H. B. 


209[L].—L. Scuwarz, ‘Zur Theorie der Beugung einer ebenen Schallwelle 
an der Kugel,” Akustische Z., v. 8, 1943, p. 91-117. 
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The pressure on the surface of a rigid sphere under the influence of a plane wave of 
sound can be calculated by the known theory of Rayleigh. Calculations were made by the 
author with the aid of Miss E. Friedrichs of the real and imaginary parts of the functions y 
and 

(2m + 1)in*t 
ye = — (2 Po w = 2xa/d = ka 

for @ = 0(5°)180° and w = 1(1)10, 5D being given. Table I gives |y| and T. II gives arc y 
for the same ranges. The first quantity represents the absolute value of the ratio |p/po| on 
the spherical surface of the amplitudes of the diffracted and incident waves, the second quan- 
tity decreased by w cos @ represents the difference in phase of these waves and is given in 
both radians and degrees. Tables are given also for the real functions fn(w), ga(w) defined by 
equation 


= (2/xw) "[fn(w) + ign(w)] 

The ranges are nm = 0(1)15 and w = 1(1)10. The tables of the quantities 

An(w) = (2m + = = (2m + 

fr = nfalw) — Bn = mgn(w) — 
are for w = 1(1)10 and for various values of ranging from 0 to 22, 6D being given. A 
table with 5D is also given for |p/po|?. Diagrams illustrate the numerical results. In 
fig. 1 relating to |p/po| it is noted that the maximum at the south pole denotes ‘optically 
speaking’ that there is a bright spot on the face of the sphere away from the source. 
Plots of |p/po| + -6w— 1, |p/po| + .040°, ¢° — -40°, ¢° — 20° against @ or ware given 
to elucidate the properties of the functions tabulated. 
H. B. 


210{L].—H. Srenzet, ‘Uber die Berechnung des Schallfeldes unmittelbar 
vor einer kreisférmigen Kolbenmembran,’’ Annalen d. Physik, s. 5, 
v. 41, 1942, p. 256-259. 14 & 21.6 cm. 


Sm(x) = (x), Cm (x) = (x) ; 
T. 1 gives the values of 
I 1-3-5- ...2m—1 é 
[Sm—1(10) + iCn—1(10)], for m = [1(1)12; 4D]; 
II. (.4)"-Jm(4), for m = [1(1)5; 4D]; 
I X II, for m = [1(1)5; 4D]; 
III. (.8)"J,.(8), for m = [1(1)12; 4D]; 
I X III, for m = [1(1)12; 4D}. 
Table 2 gives the values of 
1-3-5- ... 2m 
2-4-6- ... 2m 
II. (10/12)"*[S,,(10) + iC,.(10)], for m = [0(1)14; 4D]; 
I X II, for m = [0(1)14; 4D]; 
III. (10/16)"*[S,,(10) + 7C,,(10)], for m = [0(1)9; 4D]; 
I X III, for m = [0(1)9; 4D]. Compare SteNzEL 1, MTAC, p. 233. 


ate to 3D, for x = 10, 20, 40, 60. y = 0(1)20. For y < x 
1°3-5-...2 

+ ifm = + and for y > x 

m=1 


= 1-3-5. 2m — 


Pa + = 1 — — = + iCn—1(y)]. 


There are also graphs Of p., Pm for y = 10, 20, 40, 60; 0 < x < 20. 


+ Jm4i(10), for m = [0(1)15; 4D]; 
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MATHEMATICAL TABLES—ERRATA 


References have been made to Errata in the introductory article of this 
issue, by Bickley (Goldstein, Hidaka, Lubkin & Stoker, Morse & Ruben- 
stein), and in RMT 196 (Cayley, Kav4n), 199 (Wright), 205 (Buchholz). 


t 
62. JAMEs BurceEss, ‘“‘On the definite integral f with extended 
tables of values,” R. So. Edinburgh, Trans., v. 39, 1898, p. 257-321. 


A. The great tables of this number of the Transactions do not begin until p. 283. In 
the earlier pages are various extended preparatory values of constants. Errors in these values 
in final digits are here summarized. 

On p. 258 are 30-place values for $4/x and 2/4/x, and log (2/4/x) 


for 670, read 671 
2/s/x for 120, read 122. 


These corrections are based on values of 4/x and 1/4/x, computed to 317D and 310D, 
respectively (see MTAC, p. 200). Log (2/4/x) is entirely free from error. 

On p. 279 is a table of 31 constants and their logarithms. Of the 62 numbers comprising 
the table, each to 23D, only 12 were found to be entirely correct. Most of the errata are 
attributable to the fundamental error in Burgess’ approximation to p. In order to correct 
this table each of the values was calculated to at least 32D. 


Constants Logarithms 
for read for read 
p 351 142 30 78 28 88 
1/p 615 78 22495 69 22 7112 
a 325 1 26 61 56 5775 
pV/2 3515110381 81743 202 23 7 65 575 
2pV/4 009 806 981 30 298 912 773 92 20 15 18 25 
2p/Vx 9 82 746 8 88 698 
2p?/+/x 59 638 137 62719 00537 59598203 459 350094 49966 359 350 094 495 86 
1/(p\/x) 142 6 60 359 549 
1/(2p+/7) 071 3 30 79 85 8175 
pvr 004 903 490 65 149 456 386 96 641 451 
for read for read 

1/(p/2) 058 7 08 235 425 
prilé 301 048 599 409 
p(x /4)¥10 3 06 102 3 16 1 26 
pvV/ (x — 2) 972 750 10 25 08 35 
— 8)/6 465 247 400 210 
pv/(945e — 128)/40 50 44 47 52 957 757 
p(4/3)4 5 56 315 700 5 10 
p(4/3)"4 3 34 110 389 199 
p(113/45)4 605 274 40 23 38 33 
p(8/15)"6 61 27 59 39 90 77 88 86 
6 4 38 809 722 
592 18958800 89960812983 2362962372 8362962285 
2 06 365 191 278 
2e-P?/4/a 3 84 5 63 02 89 
e/232 42 92 

V (2/2) 15 207 88 51 207 88 


The table on p. 281 was checked by comparison of Burgess’ results with a 15-place table 
which I computed with the aid of NYMTP, Tables of Probability Functions, v. 1. The great 
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carelessness displayed in the preparation of the table is illustrated by log .08888591 given 
instead of log .088885991, and log. 272460716 instead of log .272462716. 


t log t 

H for read H for read 

A 85 991 55.990 A 832 9230 686 7124 

3 2 59 49 

4 49 59 3 3 8936 7 0795 

6 79 81 4 0986 1098 

a7 9 8 6 43 61 

9 3 4 7 85 78 
8 5 7 
9 87 90 


On p. 282 are given 30-place values of e~*, and 27-place values (with one exception) 
of 2e-*/4/x, for x = 0(1)10, and 4. I recalculated each of these values to at least 38D. The 
value of e* for x = } should end in 991 instead of 990. If for x = } the value of 2e*//x 
had been given to 27 instead of to 26 places, 367 should be substituted for 37. All other values 
in this table are correct. The value of e~* was computed to 80D and thus the value to 72D, 
given by Peters and STEIN in the Anhang to Peters’ Zehnstellige Logarithmentafel, v. 1, 
p. 12, was shown to be entirely correct. 

J. W. WreENcH Jr. 
4211 Second St., N. W. 
Washington 11, D. C. 


B. 
page t error in for read 
283 0.015 H 0.0179... 0.0169... 
0.017 53126 53113 
0.055 98209 98333 
0.110 argument 0.010 
0.115 argument 0.015 
0.155 1.0815. . 1.1015 
0.156 1.0812. . 1.1012 
0.157 /4/ 1.0809. . 1.1009 
0.158 / 1.0805. . 1.1005 
0.159 1.0802. . 1.1002 
0.160 . 50273 .. 59273 
0.188 /\/ .94388 . 94288 
285 0.291 9220558 . . .9320558 
286 0.367 H . ..0679 
287 0.4293 Ai 939383 938298 
292 Heading /4/w 
Heading A2 
0.987 Ai . 5541 . .5549 
296 1.011 Ai .615273 615373 
1.036 log [2e-/4/x] + 10 85280 . 85290 
298 1.122 log [2e-/4/x] + 10 602099 . .606210 
299 1.154 log [2e~/4/x] + 10 98049 . 98149 
301 1.250 log [2e-#/4/x] + 10 30829 . .30833 
302 1.308 ...30256 . . .30286 
1.342 Ae . . 36988 36888 
304 1.405 Ae . 406067 . ..406057 
306 1.516 log + 10 . . . . 
308 1.798 H ... 1899... 
310 1.966 H ... 7457 ... 7447 
1.998 H . ..226026 .. .326026 
315 2.492 log [2e-/4/x] + 10 7.355468... 7.355458... 
317 2.6304 Ai . . 890673 . . 490673 
318 2.782 nV) 91 1091 
321 3.5 H ... 8901628 . . 6901628 
4.1 H . . 932999724 . . .993299972 
4.7 H .980048 . . 970047 
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C. Integrand 2e~* /x+. For the following values of ¢ the last figure should be (a) increased, 
(b) decreased, by a unit: (a) .187, 1.43, 3.3; (b) .947, 1.076, 1.077, 1.112, 1.230. P. 295, ¢ = 3, 
for ...983, read ...947. 

Integral H. For the following values of ¢ the last figure should be (a) increased, (b) 
decreased by a unit: (a).397, 1.274, 1.276, 1.392, 2.504, 2.506, 2.510, 2.514, 2.552, 2.556, 
2.628, 2.630, 2.634, 2.692; (b) .886, .927, .983, 1.260, 1.347, 1.466, 2.524, 2.642, 2.658, 2.662, 
2.666, 2.668, 2.670, 2.898, 3.4. 


P. 317, t = 2.644, for ...263,read  ... 261, 

t = 2.646, for ...857,read  ... 855, 

t = 2.648, for ...153,read  ... 150, 

t = 2.650, for ...978,read ... 976, 

t = 2.652, for ...264,read ... 262, 

t = 2.654, for ...056,read ... 054, 

t = 2.656, for ...529,read ... 527, 

t = 2.660, for . 963, read . 961, 
p. 321, H, t= 6, for .. 516 075, read . _.519 737, 
p. 321, G, t=4.7, for . 544, read . 545, 

t=6.0, for 069, read 
On p. 314 the argument 2.880 should read 2.380. 

NYMTP 

Epitoriav Note: In the above lists four errors in the 96 entries on p. 321 are noted (if two 
unit errors and one two-unit error, in the last are omitted). But a report of J. W. 


WRENCH Jr. (to be published in our next issue) on the 24 entries of the L-column, shows 

that 13 are in error, several seriously. Accordingly J.O. Irwin, (BAASMTC, Mathematical 

Tables, volume VII, The Probability Integral, by W. F. SHEPPARD, Cambri a See. x), 

perhaps correctly states that this table is “Seriously infested =e “7. de whe 
given by Burgess in his footnote on p. 321, for ...483925, read . 


63. M. Kraitcuik, Recherches sur la Théorie des Nombres, v. 1, Paris, 1924. 


On p. 131-191 is Table I which gives the exponent 7 of 2 mod , for p < 300,000. For 
pb > 100,000 I have discovered the following errata (p. 159-186): 


104161 for 60 read 30 
114601 6 
121081 
127681 
267481 


A. E. WESTERN 
35 Essex St., Strand 


London, England 


64. E. C. J. v. Lommet, Bayer. Akad. d. Wissen., math. natw. Abt., Abh., 
v. 15, 1886, p. 648, T.I IIa, Maxima and minima of the Fresnel integrals ; 
also G. N. Watson, A Treatise on the Theory of Bessel Functions, 1922 and 
1944, p. 745. Compare MTE 58, p. 366f. 


The 32 values of this table have been completely checked and only the following 
four errata were found: 


S(2x/x)* C(2x/x)t 
x for read x for read 
6. 2m) .343415 .343416 10.995574(=72/2)  ..380389 .380391 
15.707963(=5r)  .600361  .600362 45.553093(=29r/2) .559088  .559087 


J. W. WrencH Jr. 
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65. NYMTP, Tables of Lagrangian Interpolation Coefficients, New York, 
1944. See MTAC, p. 314f. 


P. 391, the entry corresponding to m = 8, m = 1, and k = — 2, should be negative. 
A. N. Lowan 


66. R. M. PaGe, 14000 Gear Ratios ..., New York, The Industrial 
Press, 1942. See RMT 87, p. 21f. 


In MTE 53, p. 326f, I gave a long list of the errors in this table found by Mr. S. Jonn- 
ston. We had hoped that the list would prove to be complete, but now Mr. F. LANCASTER, 
of Huddersfield, writes that he has checked Table 4, and found the following additional 
errors: 


Page N For Read 
371 621 23 X 37 23 X 27 
388 3904 59 X 66 Delete 
391 4901 67 X 73 Delete 
393 5432 46 X 118 Delete 
400 9682 94 X 113 94 X 103 


There are also three errors of position—less serious because they are unlikely to be 


Page N 

384 2860 52 X 55 should follow 44 X 65 
398 8100 Transpose 81 X 100 and 75 x 108 
401 10192 Transpose 98 X 104 and 91 X 112 


UNPUBLISHED MATHEMATICAL TABLES 


Reference has been made to an unpublished table in RMT 202 (Bis- 
SHOPP); also to results by Ince and Bickley, MTAC, p. 412, 417. 


34[A, B].—Table of x"/n!, Manuscript prepared by, and in possession of, 
the NYMTP. 


This table is for x = 0(.05)5, m = 1(1)20, to 10S. 
A. N. Lowan 


MECHANICAL AIDS TO COMPUTATION 


15[Z].—H. P. Kueuni and H. A. Peterson, ‘‘A new differential analyzer,” 
Electrical Engineering, v. 63, May, 1944, p. 221-227. (Also in A.I.E.E., 
Trans., v. 63, 1945, and discussion p. 429-431) 20.5 X 28.6 cm. 


The article describes a differential analyzer of the Kelvin wheel-and-disc type which was 
built by the General Electric Company and put into service in Schenectady in 1943. The 
design follows closely that of the machine started in 1926 at Massachusetts Institute of 
Technology by Vannevar Bush, but incorporates a number of improvements which have 
been suggested by experience with later models, especially the one at the University of 
Pennsylvania. It has fourteen integrators, four manual input tables, and two double output 
tables; it can therefore be used for problems of considerable complexity. It is also arranged 
for operation as two independent units on simpler problems when not all of the elements 
are required. 

The most important of the design innovations is the electronic arrangement used to 
relieve the integrator disc of mechanical load, and thus to minimize slipping of the integrator 
disc with respect to the wheel upon which it rolls. The arrangement uses two beams of light 
which pass through a polaroid disc mounted upon the integrator disc and through crossed 
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polaroid discs upon the output shaft which follows the integrator. Any angular difference 
in the positions of the integrator disc and output shaft causes one of the light beams to be 
attenuated more than the other. The difference is detected by balanced phototubes and 
used to control the speed of the motor which drives the output shaft. The entire torque of the 
motor is available at the output shaft, but the torque required from the rolling disc is only 
that necessary to overcome the friction of its jewelled bearings. 

The electronic follower-system and other refinements have made it possible to operate 
the analyzer at higher speeds, thus cutting the time required for obtaining its graphical 
solutions. 


Paut L. Morton 
University of California, Berkeley 


NOTES 


37. NOTATION yer,X, yeinx.—We have been informed by J.C. P. MILLER 
that ALAN FLETCHER, of the University of Liverpool, was the inventor of 
the notation 

Y,(xi*/?) = yer,x + yei,x, 
which is used in the Liverpool Index. This note is in correction of the 
statements, MTAC p. 252, |. 23-24, and Corrigenda et Addenda p. 375. 


38. A ROOT OF y = e”.—In N 20, p. 202f. (see also N 25, p. 334f.) the solu- 
tion was given of the transcendental equation 10 log x = x or x = 10*/". 
In Assurance Mag. and J., of the Institute of Actuaries, v. 3, 1853, p. 323, 
E. J. FARREN contributes an article “‘On the form of the number whose loga- 


rithm is equal to itself.’” The equation y = e” is considered, and it is found 
that 


, 
2-3-4 2-3-4-5 


The derivation of this solution seems to have been due to ROBERT MurpPuY. 
R. C. A. 


2, 3-3 


39. Roots oF tan x = xf(x).—Compare N 18, p. 201f. In L. Cowen, 
“Alternating current cable telegraphy,”’ Franklin Institue, J., v. 195, 1923, 
p. 165f., there is a table, with 4—5S, of the first 8 roots of the equation 


tan x = 20x/(x? — 100). 
H. B. 


QUERIES 


14. TABLES OF TAN (m/n).—In some work I am now carrying on 
it is necessary to evaluate expressions of the form tan (m/n), to 10D, 
where m and » are integers, ranging from 1 to 25 inclusive. For values of 
n = 1, 2, 4, 5..., the values of the argument, and therefore the function 
may be found directly in NYMTP, Tables of Arc Tan x, Washington, D. C., 
1942. Is there a table to cover the cases m = 3,7, 11...? 

M. POoLLack 
105 Marion Road 
Oak Ridge, Tennessee 
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QUERIES—REPLIES 


15. CuBE Roots (Q 11, p. 372). In VeGa-Htisse, Sammlung mathe- 
matischer Tafeln, Leipzig, 1840, or 1849, or 1865, there is a table, p. 476-575, 
which has square roots and cube roots for x = [1(1)10000; square roots to 
12D, cube roots to 7D]. The desired function in this Query may be obtained by 
multiplying the cube roots from 1000 to 2000 by .04641 5888... = 10000-"/*, 

H. E. 
NYMTP 


Epitoriat Note: A cube-root table, of the same range as that of VEGA-HULsskE, is given in 
editions of Barlow’s Tables printed before 1930. 


16. RounpDING-orF NotaTION (Q 10, p. 335).—Devices to indicate some- 
thing of the m + 1st place in an -place table are desirable if the extra some- 
thing is occasionally useful but generally to be ignored. For this purpose the 
high and low dots of MiLNE-THomson & Comrie! (e.g., 23 < 3. < 23 
<3<RS353 <4. < 3, etc.) are to be preferred to either of the 
uses of the + sign referred to in Q 10, since the former usage leaves the last 
figure of the n-place table as it should be (i.e., rounded off). 

As an example of the possible utility of the M.-T. & C. device we may 
consider the applicability of the American Air Almanac to surface naviga- 
tion ; it is said to be used in preference to the American Nautical Almanac 
already 80% of the time, at least in the U. S. Navy. It is generally ac- 
cepted that the error of astronomical sights is probably of the order of 5 to 
15 minutes, or nautical miles, in the air, 1 or 2 only at sea. In order not to 
increase these errors the ephemerides and correction tables are given to the 
nearest minute in the Air Almanac, to the nearest tenth of a minute in the 
Nautical Almanac. It is clear in the first place that the former will satisfy 
the normal demands of sea navigation, and secondly that the greatest ac- 
curacy obtainable would be satisfied by about a third of a minute, rather 
than a tenth, in the tables. Accordingly, the introduction of the high and 
low dots into the Air Almanac,’ together with the improvement of some of 
its correction tables,* would give it the accuracy needed in the most refined 
sea navigation without affecting its convenience for air navigation. 


SAMUEL HERRICK 
University of California, Los Angeles 


1L. M. Mitne-Tuoomson & L. J. Comrie, Standard Four-Figure Mathematical Tables, 
London, Macmillan, 1931. 

2W. J. Eckert, in “Air Almanacs,”’ Sky and Telescope, v. 4, p. 12-15, 17, Nov. 1944, 
shows that the French air almanac uses such a device, but one which, in the writer’s opinion, 
is inferior to that of Milne-Thomson & Comrie, because the last figure must be altered in 
some cases, when accuracy to the nearest minute only is desired. 

3 The writer has discussed a new kind of “critical graph,”’ which is used in the same man- 
ner as a Critical table but can be read to greater accuracy, if need be, in “The Air Almanac 
refraction tables,” U.S. Naval Institute, Proc. v. 70, Sept. 1944, p. 1140-1141; Univ. Calif., 
Los Angeles, Astronomical Papers, no. 5. 


CORRIGENDA ET ADDENDUM 


P. 305, J. STEINER 1s, 1. 1, for v. 1, read v. 11. 

P. 391, 1.—16, -—15, for A = Airey, C = Comrie, M = Miller, read A = Airey, 
C = Comrie, M = Miller. 

P. 392, 1909, 28 Line 9, for x*, and =, read TI?, and I. 

P. 394, 1909, 138, 1. 2, move A.M., 18 up tol. 1. 

P. 397, 1. —9 for y, read C, and add: We here use C = In y for Euler’s constant, in place of 

the - commonly used by British writers. 

P. 403, 1. 3, for H. Béckh, read R. BOckn. 
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